Detection of branon dark matter with gamma ray telescopes 
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Branons are new degrees of freedom that appear in flexible brane-world models corresponding to 
brane fluctuations. These new fields can behave as standard weakly interacting massive particles 
(WIMPs) with a significant associated thermal relic density. We analyze the present constraints 
from their spontaneous annihilations into photons for EGRET, Fermi-LAT and MAGIC, and the 
prospects for detection in future Cherenkov telescopes. In particular, we focus on possible signals 
coming from the Galactic Center and different dwarf spheroidals, such as Draco, Sagittarius, Canis 
Major and SEGUE 1. We conclude that for those targets, present observations are below the sensi- 
tivity limits for branon detection by assuming standard dark matter distributions and no additional 
boost factors. However, future experiments such as CTA could be able to detect gamma-ray photons 
coming from the annihilation of branons with masses higher than 150 GeV. 



I. INTRODUCTION 



Dark matter (DM) is one of the most intriguing puzzles in physics. The fact that DM cannot be made of any of 
the known particles is one of the most appealing arguments for the existence of new physics. The experimental search 
for its nature needs the interplay of new collider experiments pQ and astrophysical observations. These last ones used 
to be classified in direct or indirect searches (although there are other alternatives [H [3j ) : Elastic scattering of DM 
particles from nuclei should lead directly to observable nuclear recoil signatures. On the other hand, DM might be 
detected indirectly, by observing the products of their annihilation into Standard Model (SM) particles. We will focus 
our discussion on this last alternative. 

It has been found that massive brane fluctuations (branons) are interesting candidates for DM in brane-world 
models with low tension [4] . From the point of view of the 4-dimensional effective phenomenology massive branons are 
new pseudo-scalar fields which can be understood as the pseudo-Goldstone bosons corresponding to the spontaneous 
breaking of translational invariance in the bulk space produced by the presence of the brane [5J |S] . They are prevented 
from decaying into SM particles by parity invariance on the brane. The SM-branon low-energy effective Lagrangian 
[5j |6] can be written as 

C B r = \g^d^ a d^ a - l -M 2 ^ a + ~ (4d^ a d u ir a - M 2 ir a ir a g^) (1) 

where a = 1 ... N, with N the number of branon species. 

One can see that branons interact by pairs with the SM energy-momentum tensor T^ v and that the coupling is 
suppressed by the brane tension / 4 . Limits on the model parameter from trcc-lcvcl processes in colliders are given 
by present restrictions coming from HERA, Tevatron and LEP-II, but also prospects for future colliders such as ILC, 
LHC or CLIC can be found in 6-9 . Additional bounds from astrophysics and cosmology were obtained in [J] . 

Even if branons are stable, two of them may annihilate into ordinary particles such as quarks, leptons and gauge 
bosons. Their annihilation in different astrophysical objects (galactic halo, Sun, Earth, etc.) produces cosmic rays to 
be discriminated through distinctive signatures from the background. After branon annihilation, a cascade process 
would occur and in the end, the particle species that can be potentially observed would be neutrinos, gamma rays, 
positrons and antimatter (antiprotons, antihelium, antideutcrons, etc.) that might be detectable by means of different 
experimental devices. Among them, neutrinos and gamma rays have the advantage of maintaining their original 
direction of motion. On the contrary, charged antimatter searches are hindered by the modification of the propagation 
trajectories. 
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The paper has been arranged as follows: in Section II we give a general overview for dark matter indirect detection 
by using gamma rays and we present the thermal averaged cross sections for branon scenarios. In Section III, we 
summarize the most relevant gamma-ray telescopes operating nowadays, both satellite and ground-based experiments. 
Section IV is then devoted to the calculation of the minimum detectable fluxes from different targets by using the 
already mentioned detectors. Finally, Section V contains the main conclusions about the performed analyses. 



II. GAMMA RAYS 



The differential gamma ray flux from annihilating DM particles in galactic sources can be written as [101412] : 
d$™ 1 dm 1 



-&r =^? M ^ x MiL Aa L/ l{ ' )]As (2) 

where the second term on the r.h.s. of this equation represents the astrophysical factor, i.e, the integral of the DM 
mass density profile, p(r), along the path (line of sight, l.o.s.) between the source and the gamma ray detector divided 
by the detector solid angle. On the other hand, the first term is the particle dependent part, with (<Jiv) the thermal 
averaged annihilation cross-section of two DM particles into two SM particles (labeled by the subindex i). The number 
of photons produced in each decaying channel per energy interval dN^/dE^ involves decays and/or hadronization of 
unstable products such as quarks and leptons. Due to the non-perturbative QCD effects, the analytical calculation 
of these decay chains is a hard task to be accomplished and therefore requires Monte Carlo events generators such 
as PYTHIA [13] particle physics software. By using this software and after exhaustive and statistically significant 
simulations, the gamma-ray spectra were obtained in [TT]. In these works, three different parametrizations were found 
to be able to fit all the available data in the whole accessible photon energy range for the studied channels: the first 
one for leptons and quarks (except for the top) is given by 

x i. 5 = fliexp f hxni _ hx , l2 ^^_ + ^_\ +q x i. 5 ln _ x)] x 2 -2x + 2 
dx \ x ai x cl2 J x 

the second for the top quark is encapsulated in the following expression 

and the last one for the W and Z bosons 

1.5 dA 7 „„/ l c i \ f ln bC? - x )] 



X 



The variable x stands for x = E~/M where E~ is the photon energy and M holds for the DM candidate mass. The 
concrete values for the constants in the above expressions (3)-(5) are given in [11] and are also available online in 
numerical codes [14] . The contributions of decreasing and increasing exponential factors change depending upon the 
channel. The asymptotic logarithmic terms appearing in expressions (3)- (5) are a consequence of the Weizsacker- 
Williams effect |15j . although with either additive or multiplicative behavior depending on the annihilation channel. 
For each channel some parameters are WIMP mass dependent whereas the remaining ones are constants. 

In the case of heavy branons, the main contribution to the photon flux comes from branon annihilation into ZZ 
and W + W~ (Fig. [I]). The contribution from heavy fermions, i.e., annihilation into top-antitop can be shown to be 
subdominant [TB]. In this case, the produced high-energy gamma photons could be in the range (30 GeV-10 TeV), 
detectable by Atmospheric Cherenkov Telescopes (ACTs) such as MAGIC [13 QJ] . 

On the contrary, if M < mw,z, the annihilation into W or Z bosons is kinematically forbidden and it is necessary 
to take into account the rest of channels, mainly annihilation into the heaviest possible quarks |19j as can be seen in 
Fig. [I] In this case, the photon fluxes would be in the range detectable by space-based gamma ray observatories [T5] 
such as EGRET gD] and FERMI EH H2, with better sensitivities around 30 MeV-300 GeV. 

Therefore, in order to compute the final spectrum coming from branon annihilation, we need to know the total 
annihilation cross section and its annihilation branching ratios into SM particles. For branons, the annihilation 
cross sections only depend on the spin and mass of the particle, and the expressions for branons thermal average 
annihilations (aiv) in different SM particles channels, i.e., in Dirac fermions, massive gauge field, massless gauge field 
and complex scalar field were calculated in |16j . For instance, the leading term for non-relativistic branons decaying 
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Figure 1: Branon annihilation branching ratios into SM particles. In the case of heavy branons the main contribution to the photon 
flux comes from branon annihilation into W + W" and ZZ. The contribution from heavy fermions, i.e. annihilation into top-antitop can 
be shown to be subdominant. On the contrary, if M < m\y z, the annihilation into W or Z bosons is kinematically forbidden and it is 
necessary to take into account the rest of channels, mainly annihilation into the heaviest possible quarks. 



into a Dirac fermion tp with mass m^, is given by: 



1 „ ,9 9 / „ ,9 9 \ m v 



M = Y^p M2m l ( M2 - m D\j l ~W2i ( 6 ) 



For a massive gauge field Z, of mass mz, it reads: 



(azv) = g^ys M 2 (4 M 4 - 4 M 2 m z 2 +3 m z 4 ) sj 1 - ^ ; (7) 

whereas for a massless gauge field 7, the leading order is zero: 

(<7 7 «> = 0; (8) 
and, finally, for a (complex) scalar field $ with mass m$: 



1 w2/„wa , 2\*.L TO * 2 



It is worth noting that there would be a gamma ray line from direct annihilation into photons since branons 
couple directly to them, producing a monochromatic signal at the energy equal to the branon mass. However, this 
annihilation takes place in d-wave channel and it is highly suppressed. 



III. GAMMA-RAY TELESCOPES 



For gamma-ray detection several devices have been used in the last years to observe galactic and extragalactic 
sources. Here we focus on both satellite and ground-based detectors of gamma-ray photons. 



A. Satellite experiments 



We will analyze the performance of two satellite telescopes: EGRET and Fermi-LAT. 
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1. EGRET 

EGRET (Energetic Gamma-Ray Experiment Telescope) has detected gamma rays in the interval 0.02 to 30 GeV. 
This telescope was carefully calibrated at SLAC with a quasi mono-energetic beam in the energy range of 0.02 to 10 
GeV. Although the energy range was extended up to higher energies (30 GeV) by using Monte Carlo simulations, 
we have used only data below 10 GeV [20] due to the associated larger uncertainty at high energy. The angular 
resolution was energy dependent, with a 67% confinement angle of 5.5° at 100 MeV, falling to 0.5° at 5 GeV on axis. 
The energy resolution of EGRET was 20-25% over most of its range of sensitivity. Absolute arrival times for photons 
were recorded with approximately 50 /is accuracy. The overall normalization error is usually quoted as 15 % |20j . 

2. FERMI 

The Large Area Telescope (LAT) is the principal scientific instrument on the FERMI Gamma Ray Space Telescope 
spacecraft, originally called the Gamma-Ray Large Area Space Telescope (GLAST). The LAT is an imaging high- 
energy gamma ray telescope covering the energy range from about 20 MeV to more than 300 GeV [21]. The LAT's 
field of view covers about 20% of the sky at any time, and it scans continuously, covering the whole sky every three 
hours. The LAT measures the tracks of the electron and positron that result when an incident gamma ray undergoes 
pair-conversion, preferentially in a thin, high-Z foil, and measures the energy of the subsequent electromagnetic shower 
that develops in the telescope calorimeter. The development of the reconstruction relies heavily on the Monte Carlo 
simulation of the events. The background model for this device includes cosmic rays and Earth's albedo gamma 
rays within the energy range 10 MeV to 1 GeV. Particles that might either make non astrophysical gamma rays 
and/or need to be rejected as background are included. The model does not include X-rays nor soft gamma rays that 
might cause individual detectors within the LAT to be activated [5TJ H2] • Recent investigations [53] have looked for 
potential DM subhalos targets for ACTs and FERMI at energies higher than 100 GeV claiming that multi-wavelength 
observations could play an essential role in DM indirect searches. 

B. Ground-based experiments 

On the other hand, the measurement of very high-energy gamma rays using very large ground-based Cherenkov 
telescopes is a recent addition to the variety of experiments at the interface between particle physics, astrophysics 
and cosmology. 

1. MAGIC 

In particular, the MAGIC experiment for ground-based gamma ray astronomy is the largest of the third generation 
Cherenkov telescope experiments, located at the La Palma (Canary Islands) observatory, 2200 m above sea level. It 
consists of a system of two telescopes operating in a stereoscopic mode. The stereoscopic system allows for improved 
background rejection specially at low energies (and therefore high sensitivity), improved energy and angular resolution, 
and a low energy threshold. The initial setup (MAGIC-I) consists of a single telescope (focal length of 17 m). The 
threshold for gamma detection is around 60-70 GeV with classical photomultiplier tubes (PMTs); future high-quantum 
efficiency red-extended PMTs are expected to achieve a lower threshold. MAGIC has the best light collection that 
has been attempted so far: the largest mirror with an active surface of 234 m 2 , combined with the best available 
photomultiplier tubes that can be obtained, of a quantum efficiency around 30%. As a result, MAGIC is more 
sensitive to electromagnetic showers of lower energy, and does much to close the gap existing between satellite gamma 
ray detectors (that can go up to 10 GeV energy) and Cherenkov telescopes (that presently start at energies higher 
than 100 GeV). MAGIC-I has a threshold trigger energy of ~ 50 GeV, and an analysis threshold of ~ 70 GeV at 
small zenith angle, which also permits to observe sources with higher redshift than in the past (TTJ [T5] . 

2. CTA 

We are also interested in estimating the prospects for future experiments as the Cherenkov Telescope Array (CTA) 
project, that is an initiative to build the next generation of ground-based very high energy gamma-ray instruments. 
Current systems of Cherenkov telescopes use at most four telescopes, providing best stereo-imaging of particle cascades 
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Table I: $^f^ or $^ 2 ' ) (cm 2 s 1 sr x ), where superscripts (5) and (2) denote the estimated minimum detectable flux at 5a 
or 2cr, for different detectors: EGRET, FERMI, MAGIC-I and CTA associated with Draco. The first three are current 
constraints, whereas the last one is a prospect for its sensitivity. AEff (cm 2 ) denotes the typical effective area, Att (sr) the 
angular acceptance, t exp (s) the exposure time and <E>s g (cm _2 s _1 sr _1 ) the estimated total background flux. In order to obtain 
(av) I M 2 (em's" 1 ), we have used Eq. jij with the astrophysical factor ( J}An(10 23 GeV 2 cm~ 5 sr~ 1 ) given in [25] assuming 
a NFW profile for the DM distribution. Most part of the data are taken from [25], except for the values marked with a single 
asterisk (*) which are obtained from Ref. [27] • The $b 9 signed with double asterisk (**) is calculated by means of the Eq.( 14 1 
with e = 0.01 for both MAGIC-I and CTA projects [231127] , 



over a very limited area, with most cascades viewed by only two or three telescopes. An array of many tens of telescopes 
will allow the detection of gamma-ray induced cascades over a large area on the ground, increasing the number of 
detected gamma rays dramatically, while at the same time providing a much larger number of views of each cascade. 
This results in both improved angular resolution and better suppression of cosmic-ray background events. The CTA 
will explore our Universe in depth at gamma-rays of Very High Energy (VHE), i.e., E > 10 GeV and investigate 
cosmic non-thermal processes, in close cooperation with observatories operating at other wavelength ranges of the 
electromagnetic spectrum, and those using other messengers such as cosmic rays and neutrinos. The design foresees 
a factor of 5 — 10 improvement in sensitivity in the current very high energy gamma ray domain of about 100 GeV to 
some 10 TeV, and an extension of the accessible energy range from well below 100 GeV to above 100 TeV [24] . 



IV. ANALYSIS AND RESULTS 



Although there are other possibilities, as the Galactic center, the best targets to search for a DM annihilation signal 
seems to be dwarf spheroidals (dSphs), which are the smallest known systems dominated by DM. The astrophysical 
part (J) An of the gamma-ray flux ([2| of each target depends upon the DM density. This factor is not very well 
known and introduces the most important uncertainties in these indirect detection analyses. A classic approach uses 
a Navarro-Frenk- White (NFW) profile. This profile is in good agreement with cold DM simulations and it allows an 
easy comparison with previous studies since it has been used by many authors. For this reason, we have assumed this 
profile for the Draco, Sagittarius and Canis Major dSphs and for the Galactic Center. In particular [25] : 



/OnfwM 
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r(r + r s ) 2 



(10) 



where A is the overall normalization and r s the scale radius. 

On the other hand, we have assumed an Einasto profile for SEGUE 1 since it is more consistent with observations 



PEinastoM = /9 s exp • 
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Table II: The same as Table [I] for the observation of Sagittarius. The first two columns related to EGRET and FERMI are 
current constraints whereas the last one corresponding to CTA is an estimation of its possible sensitivity. 
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Table III: The same as Table [IT] for the observation of Canis Major. 



with the scale density p s = 1.1 x 1O 8 M Kpc , the scale radius r s = 0.15 Kpc and the index n = 3.3 [T7]. Very 
recent observations [26] of SEGUE 1 (considered by many authors as possibly the most DM dominated satellite galaxy 
known in our galaxy) by MAGIC found no significant gamma-ray emission above the background when taking into 
account the spectral features of the gamma-ray spectrum of specific DM models in a supersymmetric scenario. 

In any case, both the modification of the density profile and the introduction of substructures just include an 
additional constant in the analysis that is easy to update. With our assumptions, the values of the astrophysical 
factor for each source are reported in Tables ITJV] 



We can estimate the minimum detectable flux <I> 7 taking into account the total number of observed gamma rays. 
Due to the uncertainties of these kind of analysis, it is usual to demand a significance of at least 5a. For an observed 
target presenting: exposure time of t exp seconds, instrument of effective area A e ff and angular acceptance Afi, the 
significance of the detection exceeding 5a (or alternatively 2a) is: 



cxp 
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> 5 (2). 



(12) 



The DM annihilation flux <& 7 and the background flux are given in cm _1 s _1 sr _1 [25] . The evaluation of the 
background and its value depends both on the experiment and on the source. In the case of satellite experiments, 
the diffuse gamma rays flux from astrophysical sources is the only contribution to the background depending on the 
location of the source [25] . The background for FERMI is assumed to be the same as for EGRET [25l [27l [28] : 



dNB °~ a ^GeV-W 2 - 1 - 
aE 



100 GeVV 

E B g J 



(13) 



and the exposition time t exp — 1 yr, common for the satellite experiments. We have chosen the spectral index 2.1 
since it is the most conservative value (see [2S] or [2H] for the particular values of AT for the particular targets). In 
reality, the effective area of any detector depends on the particular energy at which operates. Eq. (12 1 assumes a 
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Table IV: The same as Table [TT] for the observation in the direction of the Galactic Center. 
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Table V: The same as Table[T]for the observation of SEGUE 1. The first column corresponds to MAGIC-I and establishes the 
present constraint from this target whereas the second one estimates the CTA prospects. SEGUE 1 is the only target for which 
we have assumed an Einaisto profile for the DM distribution [IT] . 



constant value and it is the main approximation of this equation. We will assume a typical effective area for EGRET 
of 1.5 x 10 3 cm 2 , whereas we will use A e ff = 10 4 cm 2 for FERMI. On the contrary, the angular acceptance is much 
larger for EGRET: Aft = 1(T 3 sr, than for FERMI: A£l = 1(T 5 sr [DHSl 123- 

In the case of ground-based experiments, besides the above diffuse gamma-ray flux, there are two other sources of 
background: the hadronic and the cosmic-ray electrons. In any case, the hadronic source dominates at high energies 
for which the ground-based experiments are sensitive. Taking into account the data observed by the Whipple 10 m 
telescope, it is possible to find a estimation for this background rate [10l [27] : 



d* Bg -h in - 5f ,, rl -2-i -i /100GeV\ 27 

— sa e x 10 GeV cm s sr x , (14) 

AE V E Bg J 

where we have integrated above the 100 GeV threshold of MAGIC, estimating its effective area as A e ff = 5 x 10 s cm 2 
and its angular acceptance by Aft = 10~ 5 sr. The parameter e which corresponds to the fraction of hadronic shower 
which is misidentified as electromagnetic is set to the order of 1% for MAGIC [17 l H8 l 127], 

We are also interested in estimating the sensitivity of the next generation of ground-based VHE gamma-ray CTA 
instruments. Although there are still many details of the CTA project to be fixed, an important improvement in the 
effective area is expected thanks to the large number of telescopes in the array. We will assume a typical effective area 
of A e ff ~ 1 km 2 , whereas the improvement in angular acceptance and background discrimination will be typically of 



order one [23]. In Tables I|V we report the technical details of each experiment, the background estimations and the 
resulting values of the minimum detectable gamma-ray fluxes for the Draco, Sagittarius, Canis Major dSphs, for the 
Galactic Center and for SEGUE 1 respectively. 

By using the estimated minimum detectable flux at 2a or 5er significance and the particular astrophysical factor 
(J(AQ)) °f each target, the sensitivity on N\ 2 ' 5 \av) has been obtained as a function of the WIMP mass depending 
on the particular detector. The assumption of a constant effective area gives the typical power law behavior on the 
WIMP mass of this sensitivity (n\ 2 ' 5 \<jv) oc M 2 ). The corresponding curves for the different targets and detectors 




Figure 2: Sensitivity of different targets to constrain gamma rays coming from branon annihilation. The straight lines show the estimated 
exclusion limits at 5cr for satellite experiments (FERMI and EGRET). The blue thick line corresponds to the photon flux above 1 GeV 
coming from branons with the thermal abundance inside the WMAP7 1291 limits (SJcdm'' 2 = 0.1123 ± 0.0035). The area on the upper 
left corner above the corresponding lines is excluded by LEP and TEVATRON experiments for both N = 1 and N = 7, number of extra 
dimensions. 



are shown in Figs. 2 and 3. On the other hand, the theoretical value for N 1 (av) for branons has been obtained 
by integrating the differential spectrum J2i( a i v ) < ^rw L taking into account the energy threshold of 1 GeV for satellite 
experiments (Fig. 2) or 50 GeV for ACTs (Fig. 3). The resulting Nj(av) is a function of the two branon parameters 
(/, M) and it does not depend on N (number of branon species). This can be easily understood due to the fact 
that the proportionally lower flux coming from the annihilation of a larger number of branon species, is compensated 
by the higher abundance that a larger number of species provides (for a fixed coupling, i.e. for a fixed value of /). 
Assuming that the branon relic density agrees with WMAP observations [33], it is possible to obtain f(M), which 
finally allows us to plot N 1 {av) as a function of the branon mass. Thus, if the integrated spectrum line is over the 
straight lines (which represent the sensitivity at 5cr for a particular target), a detector will be sensitive to branon 
annihilation coming from a particular target. We see that present experiments (EGRET, FERMI or MAGIC) are 
unable to detect signals from branon annihilation for the targets considered. However, as shown in Fig. 3, future 
experiments such as CTA could be able to detect gamma-ray photons coming from the annihilation of branons with 
masses higher than 150 GeV for observations of the Galactic Center or above 200 GeV for Canis Major. 

It is important to note again that the above computations and figures are based on particular assumptions about the 
DM profiles and neglecting substructure contributions. Uncertainties of order one are expected for dwarf spheroidals, 
but existence of boost factors of up to three orders of magnitude has been claimed for Galactic Center analyses [30] . 
Even in this case, after having taken these uncertainties into account, only satellite experiments may have already 
observed gamma rays from branon annihilation and mainly in the low range of the spectrum. 



V. CONCLUSIONS 



We have studied the sensitivity of different gamma-ray telescopes for the observation of indirect signals of branon 
dark matter in brane-world scenarios. Under the assumption that branons are mass degenerate, this sensitivity only 
depends on two parameters of the effective theory that describes the low energy dynamics of flexible brane-worlds: 
the brane tension scale / and the branon mass M. 

We have computed the production of photons coming from branon annihilation happening in either some dSphs or 
the Galactic Center, and estimated the sensitivity for these cosmic photons to be detected in different experiments. 
In particular, we have studied the prospective detectable flux from Draco, Sagittarius, Canis Major, SEGUE 1 and 
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Figure 3: Same as Fig. [2] for ground-based detectors. In this figure, the continuous thick red line corresponds to the photon flux above 
50 GeV coming from branons with the thermal abundance inside the WMAP7 29 limits. 



for the Galactic Center for EGRET, FERMI and the future CTA. In the case of Draco and SEGUE 1, an estimation 
for the MAGIC telescope is given as well. The estimated constraints show that the interesting parameter space of the 
theory, where the thermal branon relics account for the total non-baryonic dark matter content of the Universe, has 
not been restricted by present observation yet. 

Concerning the next generation of ACTs, they seem to be able to prove this thermal area of the parameter space 
thanks to, fundamentally, the use of a large number of telescopes which can increase significantly the effective area of 
detection (Fig. [3]) . With a better sensitivity, this kind of instruments could explore the highest part of the spectrum 
for branons heavier than 200 GeV by means of the observation of both dSphs (Canis Major in particular) and the 
Galactic Center. 

Therefore, the technical instrumentation available nowadays appears to be insufficient in order to constrain the 
brane-world models through gamma-ray detection. On the other hand, the estimates for the next generation of ACTs 
show that these types of signals could provide the first evidences of these models. In the same figures ([2] and [3]), it 
is possible to see the present constraints from collider experiments. These searches are complementary and prove, in 
general, a different area of the parameter space of the model. 
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